Feeding of the elliptic flow by hard partons 
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Abstract. We propose that in nuclear collisions at the LHC the elliptic flow may 
get a contribution from leading hard and semihard partons which deposit energy and 
momentum into the hydrodynamic bulk medium. The crucial effect is that these 
partons induce wakes which interact and merge if they come together. The contribution 
to the integrated elliptic flow is estimated with the help of a toy model to about 25% 
of the observed value and shows strong event-by-event fluctuations. 

PACS numbers: 25.75.-q, 25.75.Ld 
1. Introduction 

The azimuthal anisotropy of hadron production in ultrarelativistic nuclear collisions, 
known under the term elliptic flow, has proved to be an important observable [H EJ |3] . In 
particular, its observed size lead to conclusions that concern the time of thermalisation 
[U E] and the influence of shear and bulk viscosity [HI [TJ, El El HDl EI]- To get to 
these important and interesting characteristics of the bulk matter one must, however, 
reasonably quantify the influence of other effects [121 [321 EJ. Among them there is 
our ignorance about the actual initial conditions for the hydrodynamic evolution of the 
fireball as well as the role of fluctuations in the initial conditions [TS1 HS1 Ej- Also, the 
transition from dense hydrodynamic matter to more dilute hadronic gas, which could 
slip out of equilibrium and later freezes-out, may also impact the momentum anisotropy 
[T8| H9| 120] l21j 122] . In this paper we describe a contribution to the elliptic flow which 
may be important at the LHC energy and possibly also at RHIC. 

The standard interpretation of the elliptic flow is that it is caused by the anisotropy 
of the pressure gradients within the excited matter. They are larger in the direction of 
the impact parameter, which is usually called 'in-plane' direction. (The other transverse 
direction is denoted 'out-of-plane'.) The anisotropic pressure gradients lead to different 
accelerations of the collective flow in the different directions. The fireball then finally 
expands faster in the in-plane direction. Therefore, more hadrons are emitted in this 
direction and their transverse momentum spectra are flatter than in the out-of-plane 
direction. The measured size of this anisotropy, its dependence on the transverse 
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momentum and on particle species depend on many effects. The one among them 
we want to focus at is the dependence on viscosity, since one would like to extract this 
transport coefficient from the data [T3] , 

Thus, if there is an additional cause for the azimuthal anisotropy of hadron 
production, this might influence the statements concerning fast thermalisation and 
maybe even modify the conclusion about low viscosity. In this paper, we consider 
such a mechanism. 

At the LHC, jets and minijets are produced copiously. Their energy loss when 
crossing the deconfined matter is so huge that only a few of them appear as distinguished 
jets [2U [221 EE]- They rather transfer their energy and momentum into the bulk 
matter and initiate collective phenomena there ^ZlEElESlEQlEIllSlESISlEaEB]- 
Momentum must be conserved and thus the momentum of the hard parton must be 
transformed into momentum of a stream of matter (a wake), a Mach cone wave, or 
something like this. Much interest is currently devoted to such effects. 

A question arises, what would be the result of many such streams if they all are 
initiated in the fireball? Some of these streams could merge and either cancel or flow in 
a new direction so that energy and momentum are conserved. Original hard partons are 
produced with no preferred transverse direction. The first expectation would be that 
the large number of fluid streams they initiate cancel out in some way and in the end 
there remains just thermalised matter with some energy density and no macroscopic 
flow. However, in non-central collisions the argument might not be so straightforward. 
The streams have random directions, but their spatial distribution is not isotropic, 
since it is given by the initial collision geometry. In the in-plane direction the fireball is 
narrower. Thus there is a good chance that two streams having finite width and flowing 
in the out-of-plane direction will meet (see Figure [T]). On the other hand, streams could 
more easily pass each other without interacting when produced in the in-plane direction. 
Thus — very grossly — streams perpendicular to the reaction plane cancel each other while 
those flowing in directions parallel to the reaction plane survive. As a result, the flow 
of bulk matter receives some feeding from hard partons and this feeding exhibits signs 
of a positive elliptic flow. 

Note that the scenario proposed here is related in its spirit to hydrodynamic 
simulations with hot spots in the initial conditions which are simulated on the event- 
by-event basis [151 HS1 E3 EE1 USE HQ]. Contribution to elliptic flow of the bulk from 
semi-hard partons has also been considered in (351 EE] . 

Obviously, this argument is very schematic and to make it more sound one would 
need to integrate over all possible directions of the streams. The proper way to do so 
would be the use of a hydrodynamic simulation with included energy and momentum 
deposition from hard partons. At present, such simulations are not available due to 
technical complexity. Therefore, in order to obtain simple estimates of the possible size 
of the effect, we rather construct a simple toy model, which resembles the effect. 

In our toy model which is introduced in the next Section, we represent the flowing 
streams by flying blobs of matter with specified size. If two blobs meet, they merge into a 
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Figure 1. Illustration of the likeness or unlikeness of two streams to meet. Left: two 
streams flowing in the out-of-plane direction are likely to meet. Right: as the fireball 
is elongated out of the reaction plane, two streams which flow in the in-plane direction 
have more space to pass each other without merging. 

heavier blob which moves so that energy and momentum are conserved. Originally, many 
blobs are generated. Their number and momenta correspond to the expected number 
and momenta of hard partons. In the end, after merging of all blobs which should have 
merged, they all evaporate pions. In Section[3]we then analyse the distributions of these 
pions. It turns out that about 25% of the observed integrated elliptic flow may be due to 
our effect. We also find that this contribution is very strongly fluctuating. Conclusions 
are summarised in Section HI 

2. The toy model 

Streams within the fluid are represented by blobs. For each blob first its four-velocity is 
generated. The momentum is generated according to the calculated distributions of the 
produced hard partons in transverse momentum and rapidity. For rapidity distributions 
at LHC and RHIC we assume that they are uniform in the central two units of rapidity. 

2.1. Momentum distribution 

Transverse momentum spectra have been calculated and the differential cross section 
for gluon production in proton-proton collision was parametrised as 

j^do~NN _ 1 1 daNN _ B 

dp 3 2irp t dp t dy {1+Pt/poT' 
where po, B, and n are parameters. The parametrisation works fine in the p t interval 
from 2.5 to 12 GeV/c. Calculated spectra deviate from this parametrisation for higher 
p t . Note, however, that production of jets at such high p t is rare and thus does not 
contribute much to the total yield and can be assumed to have small effect on the bulk 



when large number of collisions is analysed. For a simulation at LHC energies we chose 
B = 14.7 mbarn/GeV 2 , p = 6 GeV, and n = 9.5. 

Now we calculate the total number of blobs in a non-central symmetric collision of 
two nuclei with mass numbers A at the impact parameter b (b = \b\). The cross-section 
for the production of the leading particle with p t larger that p m is then obtained by 
integrating eq. Q 

a(pm)= [J"Z'^ dp,dy - (2) 

The mean total number of leading particles with p t > p m is then 

N .( v b) - A 2 T AA (b)a( Pm ) 

1 - (1 - T AA (b)a(p m )) 
In the last equation we introduced the overlap function 

T AA (b)= I T A (r)T A (f-b)d 2 f, (4) 

J overlap 



T A (r) =2p ^R\-r*. (5) 

where T A (x) is the nuclear thickness function. The radius of the nucleus is R A and for 
the sake of our estimates we have assumed very simple profile with constant nuclear 
density p . 



2.2. Evolution of blobs 

The blobs represent streams of bulk matter and carry the momentum of the leading 
partons. Once the momentum of the blob is given, for the simulation we need to 
determine the velocity of the blob. It will be close to c since all energy of the partons 
is basically due to momentum. Technically, we choose a very small off-shell mass of 
m = 1 GeV. Then, the velocity of a blob with transverse momentum p t and rapidity y 
is 

= {m t coshy, p t cos0, p t sin0, m t sinhy) , (6) 

where mt = ypf + m2 an d the azimuthal angle </> is generated randomly from a uniform 
distribution. 

Jets are produced in pairs which are roughly back to back in the transverse plane 
(but not longitudinally) so that the total transverse momentum vanishes. However, in 
the parametrisations of the leading parton differential cross section we have assumed 
some broadening of the transverse momentum due to intrinsic (k 2 ). Since this is an 
initial state effect, it follows that the total transverse momentum of the hard parton 
pair should be of the order (k 2 ). Following this geometry we find that the away-side 
parton is not directed precisely in the oposite way but may deviate from this direction 
by some angle a of the order 
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Here, we have assumed that (k%) -C pf. Thus the azimuthal angle of the generated away 
side leading parton shall deviate from the opposite direction by this a, which will be 
generated from a Gaussian distribution with the width yf {kf)/p^. In case of small p t 
this may become a large number. Thus if (k^)/pf > 0.7 we shall fix the width to 0.7. 

The transverse positions at which the blobs start moving are distributed according 
to the density of binary nucleon-nucleon collisions 

p b (r)=T A (r)T B (r-b) (8) 

The time at which the blob is created is delayed by some formation time r (set 
usually to 0.6 fm/c) multiplied by 7. The longitudinal position is then t v 3. Longitudinal 
rapidity is chosen from a uniform distribution. 

Two blobs merge when they approach each other in their pair CMS closer than 2Rt, 
where Rb is a model parameter which we will vary in our simulations. In this formulation 
it is the size of a blob but it actually represents the radius of the fluid stream. In the 
merger a new blob with the same size and a higher energy content is created. It will 
move with velocity v and the direction of the velocity is chosen such that energy and 
momentum are conserved. 

When there are no more mergers the blobs evaporate into pions. In the rest frame of 
the blob, pions emitted by that blob are distributed thermally with a kinetic temperature 
T. In our simulation we choose T=160 MeV. 



3. Elliptic flow 

The elliptic flow coefficient t> 2 is defined as the second order Fourier coefficient in the 
decomposition of the azimuthal single-particle distribution [2] 

P(<f>) = = N (1 + 2v 2 cos[2(0 - O )] + •••), (9) 

where <fro is the angle of the event plane. We have written this relation for midrapidity 
so that certain symmetry constraints apply. Since in our simulation we know the 
orientation of the reaction plane, we can always use the coordinate frame in which 
0o — and determine v% from 

C P(<P) cos( 20)rf0 

We measure v% in each event as the average of cos(20) over all particles and then we 
take the average over all events. 

In Figure [2] we show the azimuthal distributions summed over all simulated events. 
For the most central classes we simulated more than 150,000 events and for centralities 
over 50% the samples count 100,000 events. For the analysis we accepted pions in 
rapidity window from -1 to 1. 

In the histograms we clearly see that the production of final state hadrons is 
correlated with the reaction plane. Thus we deal with true flow signal and not just 



v 2 = Jo • (10) 
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Figure 2. Histograms of the azimuthal distributions of pions integrated over p t . 
Different panels correspond to different radii of the streams Rb and different centrality 
classes, as indicated in the panels. 



with a non-flow effect stemming from correlation of individual hadrons with each other. 
The difference between results of simulations with two chosen radii of streams is not 
very large. 

In order to better see how the anisotropy depends on centrality and the size of 
the streams we plot in Figure [3] the elliptic flow parameter i> 2 as function of centrality. 
As expected, v 2 increases when departing from central events towards more peripheral. 
This growth turns into decrease for the most peripheral collisions due to lower density 
of streams where the chance of mergers is smaller. The integrated v 2 in our simulations 
reaches up to about 0.02. Note that this amounts up to 25% of the integrated t>2 
measured in Pb+Pb collisions at 2.76 ATeV [H] and a similar share of the elliptic flow 
at RHIC (12J H3]. Thus we want to conclude that the proposed effect of momentum 
feeding from hard leading partons with subsequent stream interaction may contribute 
considerably to the observed elliptic flow. 

We also found in the simulation that i> 2 exhibits large fluctuations of the order 
of mean value or even more in some cases. This follows from the nature of the effect 
which crucially depends on the mergers of streams and vanishes if no such effect occur 
in a given event. The fluctuation might possibly be milder if complete hydrodynamic 
evolution is added to the picture. Nevertheless, we expect that one of the signatures of 
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Figure 3. Centrality dependence of v-i for different radii of the streams: i?& = 1.5 fm 
(solid red line), Rb = 2.5 fm (dotted blue line). 



the effect discussed here will be enhanced v-i fluctuations. 
4. Conclusions 

Our study qualitatively shows that in nuclear collisions at the LHC relevant effect for 
momentum anisotropy may be the transfer of momentum from the leading hard partons 
to hydrodynamic medium. Crucial ingredient is the interaction of streams induced in 
the bulk medium. They may merge and flow together. The net effect is the deflection 
of some of the streams from the out-of-plane direction into the in-plane direction. 

The toy model simulation amounted to 25% of the elliptic flow that is measured 
in Pb+Pb collisions at the LHC. The conclusion is that the effect may be important 
and should be included in simulations aimed at extraction of bulk matter properties by 
comparing to experimental data. So far, simulations with one jet depositing momentum 
into bulk matter have been performed [33J. However, to our best knowledge no 
hydrodynamic simulations exist, where the effect of energy and momentum feeding 
from many hard partons would be investigated. 
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